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TECHNICAL MEMORANDUM X-64571 


CONTRIBUTIONS TO X-RAY PHYSICS IN SPACE 

SUMMARY 


The first section of this report is a discussion of the continuous 
X-radiation in space that can originate as blackbody radiation, deceleration 
radiation, and by inverse Compton processes. 

It is shown that the radiation in Tau X-1 originates as synchrotron 
radiation, in Cyg X-2 and Sco X-1 as deceleration radiation, and the diffuse 
background radiation by inverse Compton processes. 

A statistical argument is presented to clarify the origin of the 
unidentified sources. 

In the second section the hypothesis is made that the background X-ray 
flux in the galaxy consists partly of synchrotron X rays that are generated in a 
galactic dynamo. Synchrotron radiation from the galaxy has been observed at 
radio frequencies. The radio emission of the galaxy provides evidence that 
magnetic fields extend far beyond the boundaries of the galaxy proper into a 
corona which is called the halo. The magnetic field of a cosmic dynamo is 
maintained by the induction effect of the turbulence in an electrically conduct- 
ing fluid. It is probable that turbulence exists in the halo, and therefore, the 
existence of a dynamo in the halo is also probable. 

The last section discusses maintenance of the synchrotron radiation of 
the Crab Nebula by magnetic fields. The only theory that appears to be capable 
of accounting for the required fields is the dynamo theory. The continuum of 
the Crab Nebula presents a turbulent electromagnetic field, and therefore, the 
conditions for a migratory dynamo exist. Currently, the Steenbeck and Krause 
model offers the best developed procedure to obtain magnetic fields by dynamo 
action. In this model, computer application becomes possible. It deals with 
the restricted problem, in which the velocity is considered as given. If a 
model of the Crab Nebula were known with a precision comparable to that of 
solar models, the structure of the magnetic field of the Crab Nebula could be 
computed. The radio and optical pulsar NP 0532 is also an X-ray pulsar close 
to the center of the Crab Nebula. It was shown that the principal properties of 



pulsars can be explained if the falling core of a collapsing star rebounds at 
finite or infinite density as assumed in the oscillating model of the universe. 


X-RAY ASTROPHYSICS 
Introduction 

X-ray sources in space are young objects and lie in the spiral arms of 
the galaxy. Most of them lie within several degrees of the galactic plane. 

The distribution in longitude is highly unsymmetrical. The distances for the 
X-ray sources are large, of the order of 3.084 • 10^^ m. Closer objects 
would be expected to show a larger dispersion in latitude and more symmetry 
in longitude. 

Continuous X-radiation can originate as blackbody radiation, synchrotron 
radiation, deceleration radiation, and by inverse Compton processes. 

Blackbody radiation is characterized by the Planck distribution curve 
and by very high surface emissivity. It implies sources of small physical 
area, small angular dimensions, and lifetimes too short to be acceptable for 
any great number. 

Synchrotron radiation is emitted by electrons moving in magnetic fields. 
Very high electron energies are required to produce high-energy photons by 
this process. 

Deceleration radiation or free-free emission is generated in the inter- 
action of charged particles with matter. It results from e-p Coulomb scatter- 
ing at nonrelativistic energies and from both e-p and e-e scatterings at 
relativistic energies. 

In the Compton scattering of a low-energy thermal photon by a high- 
energy electron, a high-energy scattered photon is produced with the energy 
being transferred from the electron. 

Lines cannot be expected, except possibly for iron, and will not be 
considered. 
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The following three spectral forms have to be considered: 


1. For blackbody radiation, Planck's law: 


S(W) = c 



2. For synchrotron or inverse Compton radiation, a power law: 


S(W) = CW 

3. For deceleration radiation, an exponential law: 

_W 
~ kT 

S(W) = Ce 

The problem of matching observations to such spectral functions is 
complicated by the usual occurrence that no single temperature describes a 
thermal source, and no single index fully describes a synchrotron source. 

A conclusive proof of synchrotron emission is the detection of polarization. 

The following astronomical units are used: 

1 parsec = 3.0856 X 10^® m , 

1 It yr = 9.4605 X 10^^ m , 

1 parsec = 3. 2615 It yr = 206 264. 8 AU , 

1 AU = 1.49598 X lo“ m , 

1 Ityr = 6.324 X 10^ AU , 

1 tesla = 10"^ gauss 

1 baryon = l/56 of Fe®® = 1.659 X 10“^^ kg , 

m = -2. 5 log “ . 

o 
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Generally, the International System of Units has been used. However, 
in some cases , the equations as given by the authors have been left unchanged. 

Table 1 contains the X-ray sources as identified by Fujimoto et al. {!]. 


Identified X-Ray Sources 

Tau X-l, Sco X-1, Cyg X-2, and Cen X-2 have been identified as 
X-ray sources. There is evidence that Cas A and Cas B are weak X-ray 
sources. The nature of the remaining galactic sources is not known. The 
presently observed X-ray sources are located in the spiral arms. 

The total number of X-ray sources in the galaxy has been estimated to 
about 1250. The average luminosity of the known galactic X-ray sources is 
about 6 X 10^ W. 

Tau X-l. The Crab Nebula surrounds the center of a supernova explo- 
sion which occurred in July 1054. The distance to it is 3. 639 • 10^® km. 
However, this value has been obtained on several simplifying assumptions 
(e.g. , by ignoring the ellipsoidal form of the filamentary shell). The true 
distance may be 5, 242 • 10^® km. The angular diameter of the visible 
nebulosity is 5 min. The supernova which produced the nebulosity belongs to 
type I (Table 2) . The envelope of the ellipsoidal volume is still expanding at 
about 1000 km sec“^ 

In visible light, the Crab Nebula is an amorphous mass of white gas 
within tangled filaments of red-glowing hydrogen. The white light spectrum 
is continuous and highly polarized. This polarization, in combination with the 
fact that when the X-ray data are plotted together with the data of the radio 
and optical continuum emission, the three parts of the observed spectrum 
indicate a single emission continuum — with a break extending from the radio 
to the X-ray region — is the reason for the adoption of a synchrotron model. 

Tau X-l is a powerful X-ray source and, as all known remnants of 
galactic supernova without an exception, is also a radio source. To produce 
synchrotron radio emission requires relativistic electrons with energies of 
the order of 1. 6 • 10“^“ J in a magnetic field of about lO"*^ tesla; for the 
optical emission, the electron energies have to be as high as 1. 6 • 10~® to 
1.6 • 10"® J; to produce X rays requires the electron energy spectrum to 
extend to 1.6 " 10“® J. In the X-ray range, from 1 to 10 • 10"^® m, the 
observed flux is 

L = 4 X 10^® W . 
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TABLE 1. X-RAY SOURCES [1] 
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Sources which definitely show power spectra. They are not included in calculating the number of 
X-ray sources as a function of gas mass. 



















TABLE 2. GALACTIC SUPERNOVAE [3] 


Year 

m 

max 

®1900 
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a. Rad: Radioquelle — radio source 


b. Crab Nebula 

c. Tycho's Nova 

d. Kepler's Nova 

The diameter of the X-ray source, as deduced by Friedman et al. [2] 
during a lunar occultation, is about 1. 9 • 10^® m, as compared to a diameter of 
about 5. 676 • 10^® m in visible radiation and to a diameter of about 20 km if the 
source were the hot surface of a cooling neutron star. In addition to the extent 
of the source, the shape of the spectrum shows that the observed X-radiation 
comes primarily from the filamentary part of the Crab Nebula rather than from 
a residual hot neutron star. 

The ultrarelativistic electrons emitting most of the power radiated in 
synchrotron radiation have lifetimes against synchrotron losses shorter than 
the age of the nebula. These electrons have to be resupplied continually with 
an energy rate of 10^^ W to counterbalance the synchrotron losses. Compres- 
sional waves are observed to damp out with a power loss of the same order 
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of magnitude as the power radiated in synchrotron radiation. It is supposed 
that the damping of the compressional waves acts as the energy source for the 
acceleration. Variations in the surface brightness, which are usually referred 
to as "wisps, " moving in a time scale of months are observed near the center 
of the nebula. The dimensions of a wisp may be as great as 1. 892 • 10^^ m, 
and its energy content is about 10^^ to 10^® J. These variations in optical 
synchrotron radiation are interpreted as local compressions in the magnetic 
field. The observed damping of the motions implies a loss of energy at a rate 
of about 10^^ W. The dynamo effect and the emission and reabsorption of fluid 
dynamic waves are the bases for the acceleration, which is confined to the 
central region. This damping supplies sufficient power to balance the 
synchrotron losses. In the central region, the power input to the electrons 
has to be balanced by losses caused by the escape of electrons to the outer 
region where most of the power in synchrotron radiation is radiated. 

The energy source exciting the fluid-magnetic waves and supplying 
10^^ W is centrally located and may be a thermal source in the far ultraviolet, 
although none has been observed. 

Spectral data from the Crab Nebula are fitted by the following expression: 
~ = 1. 63 • 10“^^ photons m“^ sec“^ J"^ 

3.2 • 10"^® < e < 1.76 • 10"^® J . 


For the total emission, the following values have been obtained: 
TE = 2.4X 10"l^Wm"2 (3.2 • 10"^® < e < 1.76 • 10"^^ J) , 


TE = 

1.9 X 

10"^^ W m' 

(3.2 • 

10"^® < 

e < 8. 0 • 

10"^® J) , 

TE = 

1.6 X 

10“^^ W m' 

(3.2 • 

10"^® < 

e < 8. 0 • 

10"^® J) . 


It is possible to approximate the spectrum with two power laws. 
Between 4. 81 • 10“^® and 8.01 • 10"^® J, a W~‘^’ ^ spectrum is obtained; 
between 8. 01 • 10“^® and 1. 19 • 10""^ J, a W”®’ ® variation may be fitted to the 
data. 

An X-ray pulsar in the Crab Nebula is identical to the radio and optical 
pulsar NP 0532. The X-ray pulsed flux amounts to about 5 percent of the inte- 
grated X-ray flux of the entire nebula (3 to 15 • 10*"^® m). If the total X-ray 
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flux (1 to 10 ' m) is about 3 X W then the corresponding 

X-ray pulsar flux is 1. 5 X 10~^^ W m~\ The pulsar X-ray power is about 
200 times the optical power and 2 X 10‘^ times the radio power. The X-ray and 
y-ray pulsed intensity is 7 percent of the Crab Nebula, This is equivalent to an 
average pulsed power of 1 . 6 X 10“^^ W m“^ at the top of the earth's atmosphere. 
The percentage of pulsed luminosity is of the same order as that observed in the 
soft X-ray region, indicating a softening of the pulsed spectrum between the 
optical and the hard X-ray regions. 


The principal properties of pulsars can be explained if a falling core 
rebounds at finite or infinite density, as assumed in the oscillating model of 
the universe. A possibility is that the bounce occurs in the neighborhood of 
the Planck density: 



10^^ kg m~^ . 


The intermediate configurations then become pulsating structures, and 
the rise and fall of the core is the basic regulating principle. During the col- 
lapsing phase, radiation from the envelope is comparatively weak; in the rising 
phase the core and its adjacent shells emit a blast of intense blueshifted radia- 
tion. Catastrophic collapse results in a bouncing core structure, the outer 
layers of which are supported by the periodic radiation blasts from the core. 

Cyg X-2. Cyg X-2 is a semidetached binary in which the X rays are 
generated by matter flowing from a G star through the point (the gravita- 
tional balance point in the rotating coordinate system of the binary) and being 
accreted on a compact object, which may be assumed to be a white dwarf. The 
requirement would be 10“® per year. Expected trajectories for such 

material have been computed by Kopal [4], The gas will lose orbital energy 
when it enters the region of relatively high density which should exist near the 
white dwarf. The gas will therefore lose energy and fall on the white dwarf. 
Orbital motion provides a periodic component to the observed fluxes when one 
part of the system occults another. Fluctuations in the transfer of matter from 
one star to the other can lead to irregular variations. Peaked oscillations of 
amplitude of about 0. 04 magnitude with rise and fall times of the order of 3 min, 
with intervals of approximately 11 min, are visible by inspection in optical flux 
of the candidate star. These are superimposed on a longer-term fluctuation of about 
a 1-hr period, of amplitude 0. 10 magnitude, superimposed in turn on a gradual 
decline of about 15 percent over a 5-hr interval. A computed Fourier power 
spectrum confirms these visual impressions and shows other components. The 
spectrum is definitely nonwhite. 
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There is some indication from the photometric data that the system 
may be an eclipser. The rotational motion of the object with respect to the 
galaxy must be quite small and the eccentricity of its galactic orbit quite large. 
On this evidence, Cyg X-2 must be classified as a Population II object. 

The spectrum of the candidate object for Cyg X-2 contains numerous 
sharp metallic lines identifiable with a star of spectral type G, probably a 
subdwarf. It is not certain if the short-term variations in brightness are those 
of the subdwarf or of some other source of energy. Hell in emission and Cal in 
absorption vary together in radial velocity. A short -period oscillation of 
5.7 hr is suggested by the data, but longer time-scale variations of the mean 
velocity of the oscillation suggest the presence of a third body in the system, 
although other interpretations are possible. The general behavior of Call and 
Hy suggests the presence of gaseous streams. The distance to the X-ray source 
is 1.54 • 10^® to 2.16 • 10^® km. 

In close binaries, gas streams originate when the primary overflows its 
lobe of the critical Lagrangian surface. The mass transfer will take place to the 
secondary. The energy is released largely in the form of deceleration radiation 
X rays in the 1. 6 • 10“^® J range. The total luminosity in the X-ray flux is 
about 10^ W. The only other mechanism of energy dissipation is radiation of 
gravitational waves. The gravitational radiation amounts to about 10® W if both 
components have masses of 1 M^. Variability in the X-radiation can be caused 

by changes in the rate of, or to the cessation of, mass transfer. The spectrum 
shows that the object is somewhat colder than Sco X-1 but that it possesses the 
same spectral characteristics. 

Sco X-1. Sco X-1 is starlike, smaller than 0. 5 arcsec, or 1. 496 • 10^® m, 
and is barely detectable at radio wavelengths. More than 99 percent of its total 
luminosity is concentrated in the X-ray range. 

The X-radiation is produced by free -free transitions in pure hydrogen. 

In the near-infrared spectral region, the gas is opaque and radiates as a black- 
boxy with a temperature of approximately 5 x lO^ °K, All the radiation is pro- 
duced by the same mass of gas which may be assumed to be spherical with a 
radius of O.OIR^ and a density f«10^^ electrons m“®. It requires a star with 

mass «0, 1 if the energy is supplied by gravitational infall of material. 

The best fit of the observed counts to the deceleration spectrum yields 
a temperature of 4. 8 x 10^ °K. A large part of the optical radiation cannot be 
the optical extrapolation of deceleration radiation emitted by an optically thin 
layer of a hot plasma. This is indicated by the pulsational character of the 
light variation and the well-defined color variations. The majority of the 
observations of all nights during a total time of 24 hr distributed over 1966, 


13 



1967, and 1968 forms a mean light curve with a period of 0. 5276, with an ampli- 
tude of 0. 75 magnitude in blue light, and with an asymmetry (fraction of the light 
curve occupied by the descending branch) of 0. 7. 


The distance to Sco X-1, based on the X-ray absorption, is 8.33 • 10^® m. 
This is to be considered as the minimum distance. 

The energy source is assumed to be a close binary whose streaming 
plasma transfers are accompanied by shock collisions, therefore draining 
energy essentially from the gravitational and thermal store of the stars. The 
optical emission comes from a cooler plasma volume, only partly ionized. 

The intensity at 9.61 • 10“^® J is absolute intensity (9.61 • 10“^® J) 

= 5.22 ‘ 10“^^ photons m“^ sec“^ J“^. 

The observed flux equals L (1. 6 * 10“^® - 1.6 • 10“*^ J) « 5X 10^ W, 
if the distance is about 8.33 • 10^® m. 


TE = 1. 1 X 10-1® ^ j^-2 (3^ 2 . 10"i® < e < 8. 0 • 10"l® J) . 

The observed spectrum yields a spectral flux of photons characteristic 
of deceleration radiation from a hot thin gas, thus 


1 kT 
— e 

e 

The data can then be fitted with 
kT = 6.41 • 10-1® J ^ 

The measured intensities of Sco X-1 are reproduced as given in Table 3. 

The recoil photon, which is registered, was scattered when the energy 
density of blackbody radiation was much greater than at present. However, to 
explain the diffuse X-ray background, there must have existed more or 
stronger sources of electrons in the past than at present. It is probable that 
the intergalactic gas density was greater in the past and has decreased as 
galaxies have evolved and stars formed. Such an effect would tend to increase 
the calculated X-ray background. 
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TABLE 3. INTENSITIES OE SCO X-1 X RAYS [5] 


Date 

Intensity 

(1.60-2 • 10"*^ photons 
sec“* J“b 

Field of View 
(half angles, deg) 

Accuracy 

of 

Attitude 

Calibration 

Oct 1, 1964 

4.0 

1.5 X 18.3 

Reasonable 

Pre- and post-flight 

Jun 12, 1965 

5.5 ±0.4 

10 X 45 

Poor 

Nose cone 

Jul 26, 1965 

5 ±1 

~20 elliptic 

Poor 

Preflight 

Oct 28, 1965 

3.8 ±0.3 

7 X 30 

Reasonable 

Nose cone 

Mar 8, 1966 

3.7 ±0.4 

1 X 40 

Good 

Preflight 

Jul 28, 1966 

3.7 ±0.4 

2. 5 X 8 

Good 

Nose cone 

Oct 11, 1966 

3 ±0.6 

1 X 40 

Good 

Preflight 

Feb 6, 1967 

3.3 ±0.9 

10 X 30 

Good 

Inflight 

May 18, 1967 

2.8 ±0.3 

10 X 30 

Good 

Nose cone 

Sep 29, 1967 

3.3 ±0.7 

10 X 30 

Good 

Nose cone 

Oct 2, 1967 

3.8 ±1 

5 X 30 

Good 

Nose cone 


A best fit to all the present diffuse X-ray data gives the following 
results: 

“ = 1.99 • 10 photons m~^ sec~^ J“^ sr"^ . 

The intensity of the diffuse component of cosmic X rays has been found 
to be essentially constant. The intensity at 9. 61 * 10”^® J is: 

~ (9.61 • 10“^® J) = 9.6 • 10“^® photon m~^ sec~^ J“^ sr”^ 

The energy spectrum is nearly flat between 6. 41 • 10“^® and 
1.44 • 10“^® J: 

Absolute intensity (5.77 • 10~^® - 1.44 • 10“^® J) = 3. 58 * 10^ photons 

-2 —1 -1 
m sec ^ sr ^ 
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< e 


Total emission: 7. 6 • 10 J m”^ sec“^ sr“^ for 1. 6 • 10~^® 

< 2.08 ' 10"^® J 

Total emission: 2.1 • 10”^^ J m“^ sec“^ sr“^ for 3. 2 - 10”^® < e 

< 8.0 • 10“^® J . 

There is no large change in the spectral index between 1,6 • 10”^® and 
9,61 ‘ 10“^® J. Above 1.6 • lO”^^ J, the spectrum steepens to a value of about 
2. For 1,6 • 10“^^ J < e there are indications of additional steepening. 

Therefore, the presently available experimental results for the diffuse 
X-ray flux from 1.6 • 10“^® to above 1,6 - 10~^^ J cannot be represented by a 
simple power law. The spectral index has a value of about 1. 7 from 1. 6 • 10"^® 
to 9. 61 • 10“^® J, then increases to about 2. These values are in agreement 
with the prediction of Felten and Morrison [6] based upon the inverse Compton 
effect. Their calculated values for the spectral index are 1.7 between 
1.6 • 10“^® and 1.6 • lO"^"^ J, and a larger value at higher energies. 

No dependence of the flux on galactic latitude was found within a limit 
of about 8 percent in sampling intervals of 30 deg. This result is in agreement 
with an extragalactic origin for the X-ray background. 

It is possible to explain the background of soft X-ray emission 
(A « 2 • 10~® m) by the integral emission of all the galaxies in the metagalaxy. 
The intensity measurements on Sco X-1 are consistent with the assumption that 
between 3, 2 • 10“^® to 3, 2 • 10“^® J more X-ray intensity is radiated when 
close to maximum intensity than when at minimum or average intensity. The 
change in X-ray intensity is not as great as the change in optical intensity. 

Cen X-2. Cen X-2 is within the error radius of a variable star which 
exhibits night-to-night variations with a total observed range of 0. 4 magnitude. 
The star has not undergone violent outbursts during a 60-yr optical observation. 
The distance is estimated at 1, 542 • 10^® m, which is comparable with Sco X-1, 

Cen X-2 is a highly variable X-ray star. The variability of the X-ray 
flux is shown in Table 4. The decrease in X-ray flux during the period April 
through May 1967 can be represented as an exponential decay with a time con- 
stant of 23. 4 days. This decrease was also accompanied by a softening of the 
spectrum. The Cen X-2 source gives recurring X-ray outbursts, each outburst 
lasting only for a short period of time. 
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TABLE 4. TIME VARIATION OF X-RAY INTENSITY 
FROM CEN X-2 [7] 


Experimenter 

Flight Date 

Energy in 

3.20-8.01 • 10"^® J 
Range 

(10“^^ J m“^ sec”^) 

Grader et al. [8] 

Oct 28, 1965 

<0. 25 

Harries [7 ] 

Apr 4, 1967 

11.0 ±1.0 

Cooke et al. [9] 

Apr 10, 1967 

16.0 ±1.0 

Francey 

Apr 20, 1967 

7.5 ±1.0 

Chodil et al. [10] 

May 18, 1967 

2.6 ±0.4 

Chodil et a). [11] 

Sep 28, 1967 

<0.3 

Lewin, Clark, and Smith [12] 

Oct 15, 1967 

0.62 (Extrapolated) 

Pounds 

Jun 12, 1968 

<0.1 

Rao et al. (Present Experiment) 

Nov 3, 1968 

0.68 ±0. 08 

Rao et al. (Present Experiment) 

Nov 7, 1968 

0. 83 ±0. 14 


Cas A. In the study of the motions of the luminous filaments associated 
with the radio source Cas A, the following center of expansion is found: 


Period a. (1950) 6 (1950) 

1951-1968 23. 38 deg 58. 58 deg 

Assuming a uniform rate of expansion, the explosion which generated 
Cas A is found to have occurred in 1667 AD ±8. Adopting a value 22. 5 ^ V 
for the stellar remnant of Cas A, this value is 6 magnitudes fainter than the 
value V = 16. 5 which has been obtained for the pulsar in the Crab Nebula. 

The distance of the Crab Nebula is 6. 17 • 10^® m, which may be compared with 
a distance of 1. 05 • 10^'* m for Cas A. Then the remnant of Cas A must be at 
least 5 magnitudes (100 times) fainter than the central star of the Crab Nebula. 

It may be assumed that the absorption in the direction of Cas A is very 
much greater than in the direction of the Crab Nebula. The fact that the super- 
nova outburst of 1667 AD was not observed supports the assumption that the 
absorption in the direction of Cas A is large. 
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Diffuse Background Radiation. The diffuse X-ray background has been 
established as independent of galactic coordinates to within 10 percent. The 
power law shape of the received spectrum has the form: 

1-m 
2 

m = 3.4 

between 8,0 " 10"^® and 1. 6 • 10“^^ J. The spectral shape results from the 
Compton effect where recoil photons are generated by collisions between the 
universal 3°K distribution and a leakage spectrum of cosmic-ray electrons 
from radio galaxies. However, the added flux from all galaxies within a 
Hubble distance falls short of the observations by a factor of about 10^ in 
intensity. 

The energy spectrum of the diffuse cosmic X-ray background between 
3. 2 ' 10“^® and 3. 52 • 10"^^ J fits closely a single power law: 

I (W)dW = 3. 124 • dW photons sec"^ sr"^ J“^ . 

Unidentified X-Ray Sources 

Since only a few of the known X-ray sources coincide with identifiable 
optical objects, a statistical study has been undertaken to collectively investi- 
gate the population characteristics of the galaxtic X-ray sources. The degree 
of correlation existing between the positions of X-ray sources and classical 
cepheids, galactic novae, planetary nebulae, and Wolf-Rayet stars has been 
determined. Since all these objects have low galactic latitudes, only the cor- 
relation in galactic longitude was examined; then each object forms a point on 
the galactic circle. The circle was partitioned into a number of equal intervals, 
and the correlation coefficient was calculated by comparing the number of 
X-ray sources per interval to the number of the compared type per interval. 

The results are contained in Table 5. It is seen that the correlation 
coefficient is large for the planetary nebulae and novae in comparison with 
cepheids and Wolf-Rayet stars. Assuming that, despite the uncertainty 
regarding the completeness of the surveys used, the data in Table 5 are repre- 
sentative for the actual distributions in the galactic plane of the various types 
of objects, the small correlation coefficients for the cepheids and Wolf-Rayet 
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TABLE 5. CORRELATION COEFFICIENTS [13] 


Type of Object 

Interval Size 

4 deg 

6 deg 

8 deg 

9 deg 

10 deg 

12 deg 

15 deg 

Cepheids (361) 

0.05 

0.07 

-0.01 

0.04 

0.01 

-0.03 

0.02 

Novae (92) 

0.49 

0.52 

0.63 

0. 56 

0. 55 

0.60 

0.60 

Planetaries (1024) 

0.56 

0.61 

0.69 

0.63 

0.62 

0.69 

0.67 

Wolf-Rayet (127) 

0.08 

0. 10 

0.15 

0.14 

0.32 

0.17 

0.26 


NOTE: The number in parentheses represents the objects used in 
comparison. 


stars imply that X-ray sources are not associated with these extreme 
Population I objects. The good correlation between X-ray sources and the 
galactic novae and planetary nebulae implies that X-ray sources are more 
strongly associated with Population II than with extreme Population I objects. 


Discussion of the Results 

It has been shown that the X-ray sources in space that have been identi- 
fied are either former supernovae emitting synchrotron radiation or close 
binary systems emitting deceleration radiation. The final state of a supernova 
is a superdense star. Since X-ray astronomy has advanced the knowledge of 
superdense stars and of close binary systems, the properties of these objects 
will be discussed. 

Superdense Stars. A star of 2 to 10 solar masses and density 10^® to 
10^^ kg m ^ gradually loses energy and contracts. It becomes unstable against 
collapse and finally undergoes a supernova explosion. Under appropriate con- 
ditions, the explosion leaves behind a superdense core. This star rotates, 
vibrates, and radiates. The radiations include neutrinos, X rays, long- 
wavelength electromagnetic waves, and gravitational waves. So far, no super- 
dense star has been identified. The properties of superdense stars have been 
evaluated under the following assumptions: 

1. Temperature less than 10^ °K. 

2. The matter catalyzed to the end point of thermonuclear evolution. 
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3. Magnetic induction B less than; 


2 2 
mx 


= 4 X 10^ tesla = 4 X 10^^ gauss 


0 

4. Tangential velocity of surface less than ~ . 

No continuous transition occurs from stable cold dwarf stars to stable 
cold superdense stars. The equilibrium is unstable for the range of inter- 
mediate configurations. A dwarf star can jump across this gap and be trans- 
formed into a superdense star only by a nonequilibrium process like the 
explosion of the outer part of a supernova and the implosion of its interior. 


The instability in the gap between stable dwarfs and stable superdense 
stars has its origin in the process: 


p"^ + e“ — n + r 


In consequence of this process, pressure increases less rapidly with 
compaction in stage of electron crush than it does at lower and higher pressures. 
For nonrelativistic gases over a limited range of variation, the variation of 
pressure with density can be characterized by a law of the form: 


p = C X s 

s dp 
o p ds 

A modified relation applies to relativistic gases; 


r 


o 


* + s* dp* 
p* ds* 


which reduces to the nonrelativistic result in the limit. The terms p* and s* 
are geometrical units used in general relativity. The conversion is shown in 
Table 6. 


In the Newtonian approximation, the equilibrium is stable only if 3y^ 

exceeds 4, or when the average exceeds 4/s. The gap of instability 

extends over a range of central densities from 3 X to 3 X kg m“^. 
The transition from the Newtonian theory to general relativity produces only a 
slight change in the analysis. 
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TABLE 6. CONVERSION FROM CONVENTIONAL UNITS 
TO GEOMETRICAL UNITS [14] 



Conventional 

Geometrical 

Mass of sun 

Mq - 1. 987 X 10®® kg 

M *= 1.474 km 
G 

Standard baryon 
(l/56 of Fe®®) 

At =1.659X 10"®^ kg 
s 

p * = 1. 231 X 10“®® km 
s 

Typical density 

p = 2 . 10^^ kg/m® 
e 

p - = 1. 48 X 10"^^ cm"® 
e 

Energy 

E (10"'^ J) 

M*(cm) = 0. 8258 X 10"^ E 

Time 

t (sec) 

conv 

t(cm) = ct =2. 997925X10“t 

conv conv 

Pressure 

p * 10“^newton/m® 

p* (cm"®) = (G/c^)p 


A sphere of matter of uniform density and adjustable radius r will be 
considered. A baryon number 


A = 0.74 X 10®^ , 

about three-rfpurths of the baryon content of the sun, is selected. For very 
large r, the calculated mass-energy is the mass of the corresponding number 
of Fe®® atoms, dispersed at infinite dilution, without gravitational binding. As 
r decreases, the gravitational energy and elastic forces increase in significance. 
When the density reaches 10® kg m“®, the rate of rise of the contribution of the 
elastic forces to the energy balances the rate of fall of the increasingly negative 
gravitational energy. Here, the total reaches a minimum value of -5 X 10“^® J, 
which is characteristic for the white dwarf configuration. With further compac- 
tion, the energy required for compression increases faster than the energy set 
free by gravitational attraction. At a density of 5 X lo® kg m“®, the net energy 
of interaction changes from negative to positive, and at a density of 10^'^ kg m“®, 
the net energy reaches a maximum value of 2 X lO'^'^ J, which forms an energy 
barrier. The phenomenon of the electron crush has its beginning. Consequently, 
the pressure provided by the matter does not rise as rapidly as the pull of 
gravity. At a density of 2 X lo^® kg m“®, the interaction energy turns back from 
positive-characteristic for compression to negative-characteristic for gravity. 

As the density increases, it becomes more negative. As the crushing of 
electrons onto protons to produce neutrons is almost complete, the density 
approaches nuclear values. Compressive forces regain the ability to compen- 
sate gravitational forces and balance occurs. At a density of 10^^ kg m“®, the 
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interaction energy reaches a minimum of 10“^® J, characteristic for the neutron 
star configuration at the second stable equilibrium. As the density is raised 
above 10^^ kg m~®, work done against the nuclear forces rises faster than the 
energy supplied by gravitational forces. At a density of about 10^® kg m“®, the 
net binding becomes zero. With further increase in density, the net energy 
becomes increasingly positive, similar to the energy of a spring, tightly com- 
pressed and ready to expand explosively if released. 

To summarize the result, a system of 

A = 0.74 X 10®^ baryons 

constrained to a spherically symmetric and homogeneous configuration pos- 
sesses four states of equilibrium. Two of these are stable (the uniform models 
for white dwarf and neutron stars) and two are unstable. One of the unstable 
configurations lies at the summit of a potential energy hill, separating the white 
dwarf from the neutron star configuration. The other configuration lies at the 
summit of a much higher potential hill, which separates existence as a neutron 
star from a state of complete collapse. At a certain critical baryon number, 
which is uncertain, the barrier between the neutron star configuration and the 
collapsed state disappears. For values of A of this magnitude and more, a 
stable neutron star configuration exists no longer. At another critical baryon 
number, about 1. 18 times the baryon content of the sun, the barrier between 
white dwarf and much denser configurations disappears, and no stable white 
dwarf configuration exists. 

There is no indication of the existence of a third family of stable forms 
of equilibrium in a configuration of still higher density. 

By some process of accretion, the mass of a neutron star may be 
raised above the critical limit. Achievement of supercriticality is followed by 
collapse of the neutron star which proceeds slowly at first then more rapidly. 

In gravitational collapse, the matter arrives in a finite proper time at a state 
of unlimited density. Where the matter is located and in the empty space 
around the matter for some distance, the curvature of space tends to infinity. 

A distant observer will not obtain light signals from these regions, but he will 
obtain a last ray from the surface of the infalling matter in an infinite proper 
time of his location. The surface material continues to emit radiation after 
the last ray, but these photons are caught in the constricting curvature of the 
immediately surrounding empty space. The conditions ultimately encoimtered 
by the collapsing matter and by the geometry of collapse for some distance 
outside the matter become so singular that the classical version of the general 
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theory of relativity has to be replaced by the quantura version of that theory, 
in which quantum-mechanical probability amplitudes exist for several final 
states. The range of the final states for which the probability amplitudes 
take on significant values is not known. On the time scale of the distant 
observer, the collapse does not terminate. A system in continued gravitational 
collapse shows ever-decreasing luminosity but continues to exert full gravita- 
tional attraction forever. 

A list of the galactic supernovae is contained in Table 2. 

A comparison is made between the characteristics of types I and II 
supernovae in Table?. Diagrammatic representations of the light curves of 
these types are shown in Figures 1 and 2. 


TABLE 7. TYPES OF SUPERNOVAE [14] 


Type of Supernova 

Type I 

Type II 

Ejected hydrogen? 

No 

Yes 

Predominant association with what stellar 

Old stars 

Young stars 

population? 

(Galactic halo) 

(Galactic disk) 

Presumptive mass of precursor 

1. 16 to 2M^ 

>2M^ 

Total optical energy 

4 • 10^2 J 

2 • 10^^ J 

Peak luminosity 

1036 w 

10^® W 

Typical ejected mass 

0. 1 to IM^ 

1 to 10 

Concomitant observed kinetic energy 


4 . 

Relative frequency of the two types 

less 

more 


The Crab Nebula is the remuant of the 1054 t 5 rpe I supernova. The 
pulsar NP 0532 with a period of 0.033 sec is close to its center. The period 
of the Crab pulsar is lengthening by 3. 70 X io“^ sec per day. Extrapolating 
the period back to 1054 gives 20 msec. A second pulsar has been found only 
1. 2 deg from the Crab Nebula. If the two pulsars were formed in the same 
supernova explosion, the velocity of the second must be at least 15 percent 
that of light. On the assumption that NP 0532 is a neutron star rotating as a 
rigid body, the period change implies an energy loss of 10^^ W. During this 

o -t 

period, the nebula loses about 10 W by expansion and from radiation by its 
relativistic electrons. The equivalence of these numbers suggests that the 
energy for the nebula's expansion and for the high-speed electrons comes from 
the neutron star's rotation, but it is not obvious how the rotational energy can 
be converted into particle energies. 


23 








0 25 50 75 100 IS 150 175 200 doys 

Figure 1. Schematic representation of light curve 
of type I supernova [15]. 

The X-ray flux of the pulsar NP 0532 has been found to vary by 6« 5 ±2. 7 
percent in a period of 0. 033 sec. The energy distribution of NP 0532 fits that 
of a blackbody at 10 000 °K after a correction for interstellar reddening. 


Close Binary Systems. The hypothesis is made that a binary system 
originates from the splitting of a single star into two components because of 
rotational instability during the contraction of the star. At the present state 
of the theory it is difficult to give a reliable criterion for fission rather than 
equatorial mass loss. It is known that a rotating fluid reaches a point of 
bifurcation when 


27tG s 


0. 187 . 
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It will be assumed that a star becomes rotationally unstable and splits 
into two parts when 


a = 


27 tG s 


= 0.187 


where s^ denotes the density at the center. 

The dynamical collapse of a gas cloud and the subsequent evolution of a 
spherical star until it becomes a main sequence star occurs along a Hayashi 
track. If the star rotates, its evolutionary pattern will differ from that of a 
spherical star. Once the star becomes dynamically stable, the subsequent 
contraction can be investigated assuming uniform rotation. The star contracts 
losing matter from the equator, thereby changing its angular momentum, and 
contracts down almost vertically in the H-R diagram, losing mass and angular 
momentum. During this period, it is fully convective until it reaches a point 
where the energy can be transported by radiation and a radiative core begins 
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to develop. As long as the star was completely convective, the turbulence 
distributed the angular momentum throughout the star so that uniform rotation 
was a valid approximation. With the growth of the radiative core, the inner 
regions — disregarding possible magnetic fields — are no longer coupled to the 
outer regions, and the star does not rotate uniformly even approximatelj?^ any 
longer. Each element of the radiative core contracts, therefore conserving its 
angular momentum as the star continues to lose mass from the convective 
envelope. With the uncoupling of the central region from the surface by the 
development of the radiative core, the rotational effect becomes significant 
in the central region. The parameter can then be calculated at any stage during 
the contraction, and if it becomes larger than the limit of stability 0. 187, the 
central region will become unstable and split into two parts. 

According to this theory, single and double stars are formed by similar 
processes. Then, the chemical composition in each pair was the same initially, 
but their masses may have been different. The difference in mass was respon- 
sible for the difference in the subsequent evolution. If the masses were unequal 
at the time of the formation of the system, the more massive component will 
develop an incipient hydrogen shortage and subsequent expansion ahead of its 
less massive and luminous companion. 


The Roche model of a close binary system consists of two components 
with a mass concentration so high that their gravitational potential can be 
approximated by that of a mass point. In a rotating Cartesian coordinate 
system with the origin at the center of gravity of component mj, with the 
x-axis joining the two components, and with the y-axis lying in the orbital plane, 
the surfaces of constant potential are described by the formula: 


c 


. 2 .- 1 + _2a_ 

1 + q rj 1 + q 


^2 


+ 


(1 + q)^ 


q 
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where r^ and r 2 denote the distances measured from the centers of gravity of the 
respective components. The distance between the components is taken as unit 
of length. A section of the equipotential surfaces by the orbital plane is shown 
in Figure 3, which has been drawn to scale for g = 0. 215. The two separate 
lobes coalesce in the surface containing the first Lagrangian point Lj. The 
two lobes of this particular surface represent the upper limits to the size of 
the components, if they are to be considered as separate stars. This surface 
is called the Roche limit; the two lobes are called the Roche critical lobes. 
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Figure 3. Equipotential curves of the Roche model in the orbital plane 
(Roche limit marked by heavy line) [ 16 ] . 

Close binaries are those systems in which a component fills its critical 
Roche lobe at a stage of its evolution. Even in a close binary, the components 
evolve independently of each other as long as both remain smaller than the 
Roche limit. 

A star can reach the Roche limit at phases when its outer layers expand. 
Essentially, two phases of expansion of the envelope exist; Phase A — slow 
expansion during hydrogen burning in the core and Phase B — rapid expansion 
when hydrogen is burning in the shell. 

The stages of normal stellar evolution are as follows: 

1. Premain sequence contraction. 

2. Initial main sequence. 

3. Main sequence. 
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4. Advanced main sequence. 

5. Piapid expansion. 

6. Core helium burning. 

In systems which undergo mass exchange, not all stages of normal 
evolution occur. If mass exchange occurs at Phase A, stage 4 is terminated 
abruptly, and stages 5 and 6 do not exist. 

In case of a mass exchange at Phase B, stage 5 is shorter, and stage 6 
does not occur. It is replaced by the following stages: 

1. Rapid mass transfer. 

2. Slow mass transfer. 

3. Contraction after mass transfer. 

4. The final stage of evolution of the original primary. 

The stages of large-scale mass loss are relatively short. In Phase A, 
in the system + SM^, the maximum rate of mass transfer is 

4 X Mg^y"^; in the system the maximum rate is 

2 X 10”® M^y'^ It can be expected that for Phase B, and for more massive 

stars, the rates are an order of magnitude higher. 


THE BACKGROUND X-RAY FLUX IN THE GALAXY 

Introduction 

The hjqrothesis is made that the background X-ray flux in the galaxy 
consists partly of synchrotron X rays which are generated in a galactic dynamo 
Synchrotron radiation from the galaxy has been observed at radio frequencies. 
The radio emission of the galaxy provides evidence that magnetic fields extend 
far beyond the boundaries of the galaxy proper into a corona which is called the 
halo. The halo contains at least 5 percent of the total mass of the galaxy. 
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From the rotation of the plane of polarization by the Faraday effect in 
radio sources, it follows that the magnetic fields in the galaxy have intensities 
of the order of 

gauss = 10”®^^ tesla . 

In the galactic disk, the field strength exceeds 10”® tesla; in the galactic 
halo, it is less than 10“® tesla. In the spiral arms, a value 3 • 10“® tesla is 
assumed, and in the galactic halo, a value of 6 • 10"^® tesla is assumed. In 
Table 8, order of magnitude estimates have been assembled by Woltjer [17, 18] 
for some physical quantities of the media occurring in the galaxy. A repre- 
sentative value, 10“® tesla of the induction field throughout the galaxy is given; 
this value may be too high. 

For a fully ionized gas, without a magnetic field, the coefficients of 
electrical conductivity, kinematic viscosity, and heat conductivity are given by 
the expressions: 


3 / 

CT = 6 • 10”^ T'® (Q m)“^ , 

5 / 

p = 4 • 10® n“^ T'® m® sec”^ , 

5 / 

K = 2 • 10-'^ T'2 W m“® . 


The Background X-Ray Flux 

There exists a universal X-ray background with an energy flux of about 
7.6 • 10”^^ W m”® sr“^ [19]. 

The energy spectrum of the diffuse cosmic X-ray background has been 
represented by the following power laws: 

J (e) = 3. 124 • 10“^® • e”®* ^ photons sec”^ m”® sr”^ J”^ 

3.2 • 10“‘® < e < 3.52 • J . 

The flattening of the spectrum observed at lower energies occurs below 
3. 204 • 10”^® J, with a change in spectral index of one power. A break in the 
background spectrum occurs between 1. 602 ° 10”^® and 3. 204 • 10”^® J [20]. 
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TABLE 8. PHYSICAL PARAMETERS IN THE GAS' 
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J(e) = 1.65 • 
3. 2 • 10"^® < 
J(e) = 1, 99 • 
1.6 • 10 - 1 ® < 


10”ii • f“i* photons sec”i m”^ sr“i J"! 
6 < 3.2 • 10-1® j J21] 

10“ii • e“i* ^ photons sec“i m“^ sr”i J”i 
e < 2.08 • 10“i® J [19] . 


The Dynamo Theory of Galactic Magnetic Fields 

A synchrotron origin of the background X-ray flux has been postulated 
before. Clark [22] has estimated that synchrotron X rays may constitute a 
major part of cosmic X rays incident from high galactic latitudes. He has 
calculated the expected X-ray intensity as: 

3 • 10-®^ to 3 • 10”®l J sec“l m~® sr~i Hz“i 

His estimations of the intensities of stellar and Compton X rays at 
3 • 10“1® m from high galactic latitudes are small compared to the intensity of 
synchrotron X rays. He concludes that synchrotron X rays may be an important, 
perhaps the dominant, component at high galactic latitudes at wavelengths below 
3 • 10"i® m. 


Ginzburg and Syrovatsky [23] have examined the conditions under which 
synchrotron generation may be responsible for the emission of radiation in 
the X-ray band. They discussed the case of X rays with frequencies of the 
order of v = IQI® Hz (A = 3 • 10"i® m, e = hv = 6.4 • IQ-i® J). 

This frequency is emitted only by electrons with energies 
E Rj 3. 2 • 10”® J. It is concluded that the electrons which generate this 
radiation may originate in the cosmic ray sources. 

Verma [24] has considered the possibility that the observed diffuse 
X-ray intensity in the 1. 6 • 10“*® to 6. 4 • 10-*^ J region could partly be 
explained as synchrotron emission of ultra-high -energy halo electrons. He 
estimates that the synchrotron radiation, emitted by ultra-high-energy 
electrons, in the range from 1.6 • 10“® to 1.6 • 10“* J, in the halomagnetic 
field falls in the X-ray region of about 1. 6 • 10“*®to6.4 • 10”*'* J and should 
be anisotropic. The synchrotron X rays might be distinguished from X rays 
produced by other processes, since they would be polarized if the halomagnetic 
fields were homogeneous on a large scale. The flux of these X rays would be 
larger from the direction of the galactic center than from the anticenter or the 
polar regions. Since the diffuse X-ray intensity is roughly isotropic, the 
synchrotron X-ray contribution should be small. 
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If an electrically conducting turbulent fluid is acted upon by an approxi- 
mately homogeneous magnetic field Bq , currents are generated by the induction 
effect of the turbulent motion and superimpose a perturbation field B' on the 
existing field Bq. The field B' induces in the turbulent velocity field v’ the 
electromotive forces v' X B' whose mean values are generally different from 

zero, since B' is correlated with v'. Under conditions which are generally 
satisfied in rotating turbulent systems , an electromotive force is induced in 

the direction of Bq . This induction effect of the turbulence maintains the 
magnetic field. 

The conditions in the hydrogen convection zone of the sun are known to 
some extent; this makes the computation of a solar dynamo possible. Nothing 
is known about turbulence in the galactic halo. It is probable, however, that it 
exists, and therefore, the existence of a dynamo in the halo is also probable. 


The synchrotron radiation observed at radio frequencies in the galaxy 
comes from a spheroidal halo (axial ratio 2:1) with a semimajor axis some- 
what larger than 3.700 • 10^“ m. The radiation is characterized by little con- 
centration toward the galactic center and toward the galactic plane. No stellar 
population with such a behavior is known. Therefore, it appears probable that 
the halo emission originates in a continuous field distribution and not in a large 
number of unresolved discrete sources. Irregularities in the emission may be 
due to variations in the magnetic field strength or in the density of the 
relativisitic electrons. The radiation is somewhat stronger near the galactic 
plane. This additional radiation is usually called the disk component. 

The total power emitted by one electron is given by Westfold [25]. 


P(t) 


M e^ c sin^ a 

67T ^4 



CO = 


_e_ 

m 
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For B = 10"® tesla and f = 3 • 10"^® m « 10^® sec“^, the relativistic 
parameter 

I « 6,48 • 10"® 

is required, A single electron would generate the following maximum power 
density; 

P (t) « 3.78 • 10"i® W m"2 
F (t) 3. 01 • 10"^^ W m-2 sr"i . 

The electron density in the galaxy is estimated to have the value: 

n « 4 • 10^ m"® . 

e 

Since the field is estimated to have an extension of at least 3, 0856 » 10^® m, 
the required electrons would be easily available. 


Summary 

After the known X-ray sources in space and their generating processes 
have been reviewed, a summary of the basic results obtained will be given. 

The X-radiation in the sun under nonflare conditions is generated in the 
corona as deceleration radiation of nonthermal electrons. 

The mean magnetic field of the sun is maintained by a dynamo which is 
activated by turbulence and nonuniform rotation. This field produces no 
X-radiation. A flare, however, occurs after a short-lived dynamo has been 
built up. During the first phase of the flare, the dynamo generates microwave, 
optical, and X-radiation by the synchrotron process. During the second phase 
of the flare, secondary acceleration occurs because of the collapsing dynamo 
field whose energy becomes available for this process. 

The Crab Nebula emits synchrotron radiation which is maintained by 
magnetic fields. The only theory which appears to be able to account for the 
required fields is the dynamo theory. The ex-supernova is surroimded by an 
extended atmosphere which is in a state of high-energy, turbulent motion. 

The degrees of magnetization are different for different turbulent elements. 

The continuum presents a turbulent electromagnetic field, and therefore the 
conditions for a migratory dynamo exist. 
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The following stellar X-ray sources have been identified as emitting 
deceleration (free -free) radiation: 

1. Sco X-1. 

2. Cyg X-2. 

3. Cen X-2. 


It has been determined that these three sources have the characteristics 
of close binary stars and old novae. Assuming that the explosive process has 
occurred in the nova, a shock wave will penetrate the surface forming an 
expanded atmosphere in a time which is too short for appreciable loss of radia- 
tion. When an electric field exists in a plasma with a density such that the 
electrons whose velocity surpasses a certain critical value will gain more 
energy from the field than they lose by collisions, runaway electrons will be 
generated. Runaway electrons may be the source of deceleration X-radiation in 
the atmosphere of a nova. 

No planetary nebula has been identified as an X-ray source, but a 
statistical study leads to the conclusion that X-ray sources are more strongly 
associated with disk Population II objects than with extreme Population I objects. 
Planetary nebulae emit continuous radiation which originates by free-free 
transitions involving kinetic energy losses of free electrons in the electrostatic 
fields of ions. These free -free emissions are responsible for the radiation 
observed in the radio frequency spectrum. They also contribute in the visible 
and infrared regions. Therefore, if X rays are generated by planetary nebulae, 
it can be assumed that they will be produced by free-free transitions, also. 

The diffuse background radiation results from the Compton effect where 
recoil photons are generated by collisions between the universal 3°K distribution 
and a leakage spectrum of cosmic ray electrons from radio galaxies. Super- 
imposed upon the Compton radiation is synchrotron X-radiation which is pro- 
duced in a halo dynamo. 
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FLUID MAGNETIC DYNAMO MODEL OF THE CRAB NEBULA 


Introduction 

The synchrotron radiation of the Crab Nebula is maintained by magnetic 
fields. The only theory that appears to be capable of accounting for the required 
fields is the dynamo theory. For a solution of the dynamo problem, it would be 
required to solve the fluid-dynamic and electromagnetic equations simultaneously 
to show that a motion exists which can maintain a given magnetic field and that 
this motion can be maintained by the forces available. Because of the difficulties 
of such a complete solution, the more restricted problem is usually treated to 
consider the motion as assigned and to investigate whether it can maihtain a 
magnetic field. Attempts to consider the dynamo problem in conjunction with 
the fluid-dynamic equations have been made by Taylor [26 ] and Gilman [27]. 


The Crab Nebula 

Three type I supernovae and their remnants have been recorded in our 
galaxy during historical times and are shown in Table 9 [28 ]. 

At the levels so far detectable, only the Crab Nebula is an X-ray source. 
The ex-supernova, a collapsed star, is surrounded by an extended atmosphere 
which is in a state of high-energy turbulent motion, the degrees of magnetiza- 
tion being different for different turbulent elements. An average of W is 
injected from the nucleus into the Crab Nebula, which corresponds to «10^ 
relativistic electrons per second. This quantity has to be produced in the 
nucleus of the nebula. 

The Crab Nebula consists of two different parts; a net of thin gaseous 
filaments arranged in the form of a shell at the periphery of the nebula and a 
continuous emission which occupies the entire volume of the nebula. 

The gaseous shell ejected in the supernova explosion is identical with 
the system of filaments. The filaments represent sharply defined features, 
whose emission is concentrated in individual lines. Seen in Ha, the shell is 
3.784 • 10^® m wide and 5, 68 » 10^® m long. The principal physical parameters 
describing the conditions in the filaments are: 
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TABLE 9. THREE TYPE I SUPERNOVAE AND THEIR REMNANTS [28] 



Supernovae 


Name 

Crab 

Tycho's 

Kepler's 

Year 

1054 

1572 

1604 

1° 

184.6 

120.1 

4.5 


-5.8 

+1.4 

+6.8 

m (max) 

V 

[-5?5] 

-4’^0 

-2?2 

Interstellar absorption 

1^6 

3^3 

[2“"] 

Distance 




(3. 0856 . 10*® m) 

1500 

2400 

9900 

log radio flux at 
100 MHz in units of 




10"2fiWm“® (Hz)"* 

3. 25 

2.36 

1.90 

Spectral index 

-0.26 

-0.61 

-0.61 

Radio diameters 

3'. 5 (E.-W. 
nearly Gaussian) 

6' (ring 
source) 

3'.0 (small- 
scale structure) 

log absolute radio 
luminosity at 
100 MHz [10-^ W 




sec~* (Hz)“*J 

24.68 

24.20 

24. 97 

Expansion velocity 




(km/sec) 

1400 

5900 

12 100 


R « 2.5 • 10^® m 

N « 10® m-® 
e 

M « 0. 1 
T « 17 000 °K 

B w 3 • 10"^ gauss = 3 * 10”® tesla . 
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These conditions do not apply to the region where the synchrotron 
radiation is generated. In the continuum, the nebula possesses structure which 
has no relation to the filaments. The continuum is resolved into fibers. An 
explanation for the existence and alignment of the fibers is to be found in the 
confinement of a current filament by its own magnetic field, which is known 
as the "pinch effect. " The cause of the constriction is the electromagnetic 
attraction between parallel currents. This constriction of currents is 
accompanied by an accumulation of matter, which explains the fibrous 
structures. 

The filamentary shell encloses a region emitting radiation in the optical 
continuum, which displays varying degrees of polarization with position. Com- 
pared to the mass of the shell, the continuum is negligible, of the order of 
10“^ Mq. The energy spectrum can be expressed as a power law: 


d e e 
/3 = 1 -.2 Q! . 

The term a denotes the spectral index, which changes near, 10^“^ Hz from about 
-0.3 at low frequencies to -1.0 at high frequencies. The electron-energy 
spectrum is therefore roughly from e » 0 to e « the energy of 
particles emitting near 10^^ Hz, and E"^ from to some upper cutoff energy Cj. 

The X- and y-ray spectra can be joined smoothly to the optical data. 

The ultraviolet region of the spectrum has not been observed, but the correct 
ionization conditions are produced in the filaments if the ultraviolet synchro- 
tron radiation is taken to be an interpolation between the optical and X-ray 
regions according to the power law. Therefore, the radiation from the Crab 
Nebula can be considered as a unified synchrotron spectrum that extends from 
about 10® to 10®° Hz. Continuous injection of a power-law electron spectrum 
at a rate proportional to (t - t)°’®°, where t denotes the age of the nebula 
approximately, results in a spectrum of the type observed. 

Near the center of the nebula are local enhancements of the general 
synchrotron emission, which are unrelated to the filaments. The features 
undergo rapid changes. A characteristic configuration, however, recurs and 
is called a wisp. The wisps emit a continuum which is strongly polarized. 
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The picture of the continuum can be summarized as follows: 

1. The continuum shares the general expansion, demonstrated most 
clearly by the filaments. 

2. Superimposed on this expansion is a field of chaotic motions 
throughout the nebula, but with a trend toward decreasing velocity with increas- 
ing distance from the center. 

3. The fastest motions are confined to the central area. 

The wisps are also local compressions in the magnetic field. The 
direction of propagation is perpendicular to the local field direction. Near the 
center of the nebula, the wisps propagate at about 0,6 * 10® m, damping out 
after five or six periods. Further from the center, the enhancements propagate 
at about 0.06 • 10® m. The observed translations of the brightened regions are 
formed by mass motion, in which the gas and the field move together. The 
motion is quasi -harmonic, the restoring force being supplied by the magnetic 
field. The thermal particle density of the main wisp is about: 

» 10® m~® . 

th 

This number is six orders of magnitude less than the known densities in 
the filaments. The physical picture is that wisp 1 plays the role of a piston 
moving back and forth, generating compressional waves which propagate 
outward. 

In the central region of the Crab Nebula, the field lines are generally 
aligned and parallel to the wisps, although the wisps sometimes curve outward. 
This is the configuration which might be produced by a small object emitting 
ionized gas in all directions in an initially uniform magnetic field. The particles 
emitted perpendicular to the initial field drag the field and plasma outward, 
leaving behind a hole devoid of field and the original plasma. The interface 
between the hole and surrounding plasma forms a compressed region, which is 
identified as wisp 1. Motions result from lack of perfect balance between the 
outward particle pressure and the inward magnetic pressure. The motions in 
turn generate compressional waves travelling outward from the center. This 
model is basically symmetrical, but some departure exists from exact 
symmetry. 



We are looking at a surface of rotation perpendicular to the axis. The 
power output is presently about; 

L « 7 • 10®° W . 

The total energy injected into the nebula in the form of relativistic 
electrons is: 

W « 5 ■ 10^® J . 

The average magnetic field is: 

B « 3 • 10"® tesla . 

Near the wisps, the Alfven velocity equals; 

V. » 3 * 10^ m . 

A 

The ex-supernova is surrounded by an extended atmosphere which is in 
a state of high-energy turbulent motion. The degrees of magnetization are 
different for different turbulent elements. The magnetic field of the nebula is 
essentially radial between the renanant and the expanding shell. The continuum 
presents a turbulent electromagnetic field, and therefore, the conditions exist 
for a migratory dynamo. 

The particles Injected into the field are preferentially relativistic 
electrons. The contribution from relativistic protons and heavy nuclei is 
insignificant. 

The radio and optical pulsar NP 0532 is also an X-ray pulsar. It is 
close to the center ofjhe Crab Nebula and was formed in the supernova 
explosion of 1054. 


Magnetofluid Dynamos 

A number of promising dynamo theories have been proposed, almost 
always with a view to solar application. None of them can be applied to the 
Crab Nebula numerically yet, since the structure of the Crab Nebula is not 
known with sufficient accuracy. A survey of the available stellar dynamo 
models will l)e given, stressing the basic ideas, not the mathematical 
application. 
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Parker's Model. Parker's migratory dynamo is intended to explain the 
sunspots and some of the solar prominence activity f29]. The generation of a 

magnetic field B by a fluid with velocity v and conductivity a is described 
by the induction equations: 


II 

CDl 

curi 

(v X 

B) + ^ 

fx a 
o 

V^B 

II 

(v X 

B) 

curi 

fx a 

B . 


o 


By a suitable Cartesian coordinate system and suitable boundary condi- 
tions, equation (1) is reduced to the form: 


8t 


9A 

9t 


TT 9 A 1 

H — + V^B 

9 X p CT 
o 


PB + — V^A 
IX a 
o 


( 2 ) 


The term P has the dimension of a velocity and is a measure of the 
violence of cyclonic fluid motions; P and H are assumed to be constants; and 
equation (1) is reduced to two simultaneous linear equations [eq, (2) ] with con- 
stant coefficients. Equation (2) is solved by the migratory dynamo waves: 


B = 


_ 1 (wt + Kx) 

B e 
o 


A= 

o 


(3) 


In equation (3) co can be expressed in terms of the constants in 
equation (2). 

Parker uses these waves as explanation for the fact that sunspots appear 
in low latitudes only, for the migration of the region of formation of sunspots 
toward the equator, and for the reversal of spot polarity between cycles. 
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Braginsky's Model [30 J. If the magnetic Reynolds number of the fluid 
is large, 

1 « R = n a VJl , 
m o 

a very weak generation process is sufficient to maintain the magnetic field. 
According to Cowling's theorem, no generation process exists in the case of 
axial symmetry; therefore, if 1 « R^, it is of interest to consider the case 

of a small deviation from axial symmetry, where a generation process may 
exist which, although weak, may be sufficient for dynamo operation. At first, 
the induction equations are expressed for the case where the velocity and 
magnetic fields have axial symmetry. To separate the axially-symmetrical 
part of any quantity from the purely variable part, the operation of averaging 
over the cylindrical coordinate <p is introduced. An equation for the part of 

B variable in (p is obtained by subtracting the inductions averaged over cp 
from the original induction equations. The attempt is made to obtain a solution 
of the resulting equation in the form of an expression in powers of R The 

^ 9 B 

process of field generation is then represented explicitly in the terms for — and 
9 A 

— . In the case of inhomogeneous rotation, the second terms can be inter- 
9t 

preted as the drawing out of the magnetic field lines from the meridional field 
in the azimuthal direction, and the azimuthal electromagnetic force generates 
a meridional from the azimuthal field. 

Taylor's Model. Although Taylor applied the self-excited dynamo to the 
earth, it may also be applicable to the Crab Nebula. Taylor's model uses, in 
addition to the induction equation (1), the equations of fluid dynamics in the 
form; 


2p (oj X v) = (J X B) - Vp 


div V ^ 0 


(4) 


|f = -(Tv)p . 

To solve these equations, Taylor proves the following theorem; If the 
rotating fluid is contained in a rigid envelope in the form of a surface of 

revolution, then a velocity v compatible with the equations (1) and (4) will 
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exist if and only if the magnetic field satisfies the following constraint, 
expressed in cylindrical coordinates (r, cp, z): 


/ 

r=const 


(J X B) d(p dz = 0 
(P 


(5) 


To show that equation (5) is necessary, the (p-component of the first 
part of equation (4) is integrated over that part of the cylinder which lies 
within the fluid, considered to be enclosed in a sphere, but the proof applies to 
any surface of revolution. The integrand on the left-hand side is expanded, 
and finally transformed according to Gauss's theorem. The necessity of the 
condition follows then from the second part of equation (4). 


The proof that equation (5) is sufficient is given by providing an 

algorithm for the determination of v in terms of B and p . If equation (5) 
is to be sufficient, it has to be shown that the first two parts of equation (4), 
together with the boundary condition, 

(T . = 0 , 

on the surface of the sphere, determine v if B and p are given. 


Equation (5) expresses the fact that the moment exerted by the magnetic 
forces on any cylinder of fluid, coaxial with the axis of rotation must vanish. 

If equation (5) is satisfied, the instantaneous fluid velocity is specified by the 

instantaneous values of B and p , together with the boundary conditions of 
vanishing normal velocity over a sphere, or other surface of revolution satisfy- 
ing the conditions of the theorem, except for an azimuthal velocity u (r) . 

This azimuthal velocity is determined by the requirement that equation (5) con- 
tinues to be satisfied as B develops according to the electromagnetic 
equation (1). 

Gilman's Model. A Rossby wave is defined to be a nearly, but not 
entirely, horizontal wavy or eddying flow pattern in a rotating fluid in which 
Coriolis forces nearly balance horizontal pressure forces. It is assumed that 
a latitudinal gradient large enough to excite Rossby waves exists in the con- 
vection zone of the sun, and the fluid magnetic effects of the waves are examined. 
It is also assumed that the Rossby number is small: 


R 

o 


V 


2Ql 


« 


1 . 
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The basic equations of naagnetofluid dynamics are introduced in dimen- 
sionless variables. The equations are the horizontal and vertical equations of 
motion, the equation of continuity, the horizontal and vertical magnetic field 
induction equations, the magnetic field continuity equation, and the heat con- 
duction equation. The dependent variables (pressure, density, specific volume, 
and temperature) expressed in power series in terms of are assumed to 

be small. If K denotes any of the variables, it is expanded in the form: 


K = K 


( 0 ) 


+ 


R K 
o 


( 1 ) 


R 

o 


+ 0 (R ) 


and substituted into the equations of fluid-dynamics. Collecting first-order 
terms, it is noticed from the horizontal magnetic field induction equation, 
which is the prediction equation for the toroidal field, that it contains two terms 
representing the stretching of poloidal into toroidal fields. From the first- 
order terms of the vertical field induction equations, and including a second- 
order term, the induction of the poloidal field from the toroidal field — 
accomplished by the vertical motion twisting toroidal field lines up into the 
vertical — ■ is obtained. 

Four first-order and one second-order equations form the equations for 
the dynamo. As in Taylor's dynamo, the motion is not given but is generated by 
latitudinally nonuniform heating. The induced magnetic fields react upon the 
inducing motions. Therefore, a state of balance should be reached in which, 
on the average, induction of new fields is balanced by dissipation. 

The Steenbeck and Krause Model. Currently, the Steenbeck and Krause 
model offers the best developed procedure to obtain magnetic fields [31]. In 
this model, Parker's [29] and Braginsky's [30] ideas are utilized and expanded 
such that computer application becomes possible. It is concerned with the 

restricted problem in which v is considered as given. If a model of the Crab 
Nebula were known, with a precision comparable to that of solar models, 
the structure of the magnetic field of the Crab Nebula could be computed. 

The first part of equation (1) and the second part of equation (2) are 
used and expressed in spherical coordinates. The terms A and B are developed 
into a series of Legendre polynomials. By considering the relations existing 
between Legendre functions, an infinite system of ordinary differential equations 
is obtained. This system decomposes into two uncoupled systems, and the 
problem to find nontrivial solutions for the two systems that are regular at the 
origin and satisfy the boundary conditions at the surface leads to a two-parameter 
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characteristic value problem for two characteristic values. This problem for an 
infinite system of ordinary differential equations is solved in such a way that it 
is broken off at a particular value N. The remaining finite system is solved. 

To test the convergence, N is increased step by step until the characteristic 
values are unchanged within the required accuracy. The integration of the 
system of differential equations is performed in accordance with the Runge-Kutta 
procedure. The solution of the characteristic value problem is performed with 
an electronic computer. 


The Bouncing Core Theory of the Crab's Pulsar 

It was shown by Harrison [32] that the principal properties of pulsars 
can be explained if the falling core of a collapsing star rebounds at finite or 
infinite density as assumed in the oscillating model of the universe. The 
condition that the metric remains nonsingular — with the exception of the 
infinite density bounce — can be satisfied. The intermediate configurations 
then become pulsating structures, for which the rise and fall of the core is the 
regulating principle. In the rising phase, the core and its adjacent shells 
emit intense blueshifted radiation. The periodic emission and the accurate 
timing of pulsars are explained by the theory that catastrophic collapse results 
in a bouncing-core structure, the outer layers of which are supported by 
periodic radiation from the core. 

The observed pulsed emissions originate from the underlying layers, 
which participate in the motion of the core but oscillate with limited amplitude. 
The pulse of enhanced radiation emitted during the rising phase is seen by a 
distant observer to be contracted in time. Between pulses, the radiation is 
reduced and the time is dilated. This is because the observed pulse is short 
compared with the period. Radiation from underlying shells adjacent to the 
core will also augment the pulse, but this radiation will cease between pulses. 

Energy is radiated at the expense of increasing the binding energy of 
the original core mass at infinite dispersion. In the bouncing-core theory of 
pulsars, the period is regulated by gravitational collapse and expansion. The 
observed period is increased by redshift. The central region has no equilibrium 
configuration. The pulse originates as a burst of radiation from shells adjacent 
to the core during expansion. Its observed width is contracted because of the 
predominating Doppler effect. Pulsed emission is an essential characteristic 
of bouncing-core configurations. The energy is gravitational in origin. Loss 
of energy causes the ceiling radius of the core to shrink towards the 
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Schwarzschild radius and increases the period. The pulsars observed so far 
are approaching the end point of stellar evolution and are expending the last 
remnants of available energy. 

For NP 0532 in the Crab, for IM^ , the rate at which energy escapes 

(y 

according to the theory is 6 • 10^^ W , in rough agreement with the energy 
requirements of the nebula. 

Most supernovae are possibly poorly-clad bouncing cores and are 
therefore short-lived pulsars relasing energy rapidly. 

To obtain nonstationary solutions of general relativity, the following 
assumptions are made: 

1. The field equations are valid in their most general form, containing 
the cosmological constant; 

«ik ® ^ *ik “ '"ik • 

2. The velocity of matter is negligible in comparison with the velocity 
of light. Then the momentum-energy tensor is given by: 

T., =0, if i, k 4 
ik 

T 44 = C^ p g 44 . 

3. The metric is given in the form: 

22 22 

ds^ = (dx^ + sin^ dx^ + sin^ x^ sin^ x^ dx® ) + dx^ 

with R and M depending on x^ . By a suitable choice of x^ , the condition 

M = 1 


can be satisfied without loss of generality. The metric can be brought into the 
form: 


dr' 


ds^ 


c2 


dr ^ = - (dx^ + sin^ x^ dx^ + sin^ x^ sin^ x^ dx^ ) + dx^ 
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The terms R and p are to be determined from the field equations. 
Fori = k = 1, 2, 3, equation (6) yields a I’elation; 


2 

2RR^^ cf_ 
R2 ^ R2 R2 


- A = 0 


( 7 ) 


The integration of equation (7 ) , writing t for x'^ , yields the following 
equation: 



A - R + 


A 

3c2 


R 




( 8 ) 


The integration constant is denoted by 
obtained by inversion of an elliptical integral, 


A . From this equation, R is 
i. e. , by solving the equation: 


t = 



dx + B 


(9) 


with respect to R. Constants are denoted by B and A , For the nonstationary 
solutions, equations (8) and (9) have to be considered. The term A remains 
undetermined; it is assumed that A can take arbitrary values but lies in the 

, In that case, R becomes a periodic function of t , 

with the period t * t given by the formula: 

7T 7T 


interval 1- ” , 


4c‘ 

9A‘ 



The term Xq is a function of A and A : 
X(, = © (A,A) . 


The radius of curvature varies between 0 and Xg . This is called 
the periodic solution. For small A , the period is represented by the approxi- 
mate formula: 


t - 

7T 


7tA 

c 
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